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Abstract 



We examine the effects of a van Hove singularity (vHs) in the density of 
states (DOS) of a two dimensional, t — t' — U Hubbard model within the local 
approximation. Within our approximation the non-interacting Fermi surface 
is always retained away from half-filling; and we find that the vHs in the 
DOS survives, and for hole doping is nearly "pinned" to the Fermi level due 
to the enhancement of the effective mass, m* . We discuss the electron-hole 
asymmetry of this phenomena and its relevance to the electron and hole doped 
cuprates. At optimal (hole) doping, with the Fermi level at the vHs, we find 
that the single particle damping rate F, and the resistivity pdc, scale as T 
for moderately low temperatures (T > lOOK), while for the electron doped 
system F and pdc ~ T'^- 
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The cuprate superconductors are generally regarded as doped "Mott" insulators. The 
parent state is an antiferromagnetic insulator, with a large local moment and a gap ~ 0(2ev) 
and cannot be understood within a one particle, band-theoretic approach. Transport 
measurements on the hole doped systems indicate a "carrier density" which scales with the 
doping 6. Naively one might take this to imply small hole-like Fermi pockets with area ~ 6. 
Angle resolved photoemission (ARPES) however reveals 0] a large Fermi Surface (FS) with 
area ~ 1 — 5 in agreement with the Luttinger theorem, and a shape roughly consistent with 
LDA calculations. Thus, the one particle spectrum in the hole doped system has band-like 
features, while the transport is "anomalous" in its doping and temperature (T) dependence. 
On the other hand there are systematic differences between the hole doped {p type) and 
electron doped {n type) systems. The p type cuprates exhibit dc resistivity, p^c oc T at 
optimal doping (maximum Tc), have Tc ~ (40 — 150)K and a strongly anisotropic gap 0] 
in the superconducting phase. The n type cuprate Nd2--5Ce5Cu04_j/ (NCCO) on the other 
hand shows pdc oc T^, has a rather low optimal ~ 257^^ and an isotropic gap 0. Analysis 
of ARPES data provides an interesting clue 0] to what might well be responsible for this 
difference between p and n type cuprates. ARPES clearly reveals the presence of flat 
bands near (^, 0), i.e a saddle point in the "dispersion", within 10-20 meV of the Fermi level 
(Ep) for p type cuprates. The saddle point is seen in NCCO also, but ~ 350 meV below 
Ep, and thus cannot play a significant role in transport. This suggests a possible connection 
between the anomalous T dependence of pdc and the proximity of the saddle point to Ep. 

The presence of a logarithmic van Hove singularity (vHs) in the density of states (DOS), 
arising out of a saddle point in the quasi 2-d dispersion, had been invoked earlier |^ as a 
possible mechanism for the "high Tc" and also as an explanation for the anomalous transport, 
but within a weak coupling approach [0] . While the presence of a vHs in the bandstructure of 
a quasi 2-d system is understandable within LDA calculations, the survival of this feature in a 
strongly correlated system close to the Mott transition is not, and has not been demonstrated 
in any controlled theory so far. We address the problem of the interplay of bandstructure 
and strong correlation effects in this paper, by studying a. t — t' — U Hubbard model in 
d = 2, with a vHs in its DOS (at ~ 0.2 hole doping) and with U sufficiently large to drive 
the half-filled {6 = 0) system insulating, using the local approximation. 

Within our approximation for the electron self-energy (detailed below) which retains the 
non-interacting FS at T = 0, 5 7^ 0, we find the following results (i) The vHs survives in 
the DOS of the interacting system, and its shift with respect to Ep with doping is reduced 
roughly as AE{6)/m*{S) (Fig. 1), where AE{6) is the shift in the noninteracting system, 
and m*{6) is the effective mass enhancement due to correlation, (ii) The "dispersion" of the 
interacting spectra (Fig. 2) is reduced by a factor ~ 3 with respect to the noninteracting 
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t — t' band dispersion. The resulting "flat-bands" are close to Ep for the optimally doped 
p type cuprates, and much farther down for the n type cuprates. {Hi) The single particle 
damping rate F, (Fig. 3) and pdc (Fig. 4) scale at moderately low temperatures as T for 
optimal hole doping, while they scale as for the electron doped case, {iv) However, this 
approximation fails to reproduce the temperature dependence of the Hall coefficient, Rh-, 
although the 5 dependence is recovered. 
We use the Hamiltonian 



where = —2t{cos{kxa) + cos{kya)) — At' {cos{kj.a)cos{kya) — 1) is the dispersion of the 2-d 
tight binding square lattice with nearest neighbour hopping t and next nearest neighbour 
hopping t' (the shift 4t' maintains the band edges at ±4t). We set the half-bandwidth 
D = 1, t = -D/4, t' = — 0.3t so that the vHs is located at the Fermi level for ~ 20% hole 
doping, and further choose U/D = 4. The t — t' model provides the simplified "universal" 
2-d bandstructure of the cuprates, while the choice U/D = 4 ensures that the half-filled 
state is insulating; the t — t' — U model can thus be regarded as the minimal model for the 
cuprates. 

There is no controlled method for studying the doped phase of the strong correlation 
problem defined above. We adopt an approach which is exact for the t — U model in the 
limit when the dimensionality (i — oo |^, but is necessarily an approximation for the 2-d 
problem we want to study; we assume that the self-energy and all the irreducible vertex 
functions are local, i.ek independent. In this approximation the Hubbard model maps on 
to a self consistently embedded Anderson impurity problem . A recent method devised 
by Kajueter and Kotliar [|l^] allows an interpolative solution of this self-consistent impurity 



problem for arbitrary U and 5 at T = 0, which is exact in the atomic and band limits, and 
makes feasible the study of the real frequency single particle spectra. An important feature 
of this solution is the imposition of the Friedel sum rule for the impurity problem, which 
translates to the Luttinger sum rule for the corresponding lattice problem, ensuring the 
correct low energy structure in the spectra, as verified [1^ by comparing with exact results. 
In an earlier paper |Tl[] we have used a natural extension of this method to finite T, in the 
context of a symmetric, semicircular DOS for the non-interacting band, to study transport, 
in particular the Hall coefficient, in a doped Mott insulator. We refer the reader to these two 
papers ll^jinil for the details regarding our approximation scheme. The crucial difference 
between [0] and the present study is in the consistency relation for the local propagator, 

= J depoie)/iu + p-e- S(^)). (2) 
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We now use a po(e) which closely mimics the features of the DOS of the 2-d t — t' band, 
namely the analytic form po(e) = ANln((e — ^vhY + if), with |e| < D, where = — 0.3-D 
is the location of the vHs (as in the actual t — t' DOS), t] = 10^^ D regulates the divergence, 
and An normalises the DOS. The regularisation makes the numerics simpler, and in actual 
systems can be thought of as arising from weak interplanar coupling or impurity scattering 
effects. 

We have studied the self-consistent problem for five representative choices of S; 0.08 (un- 
derdoped), 0.20 ("optimally" doped), 0.40 (overdoped), and - 0.20 and - 0.40 (electron 
doped), for T ~ O.OID — O.ID. Since the irreducible vertex functions are assumed to be k 
independent, vertex corrections to the conductivity are absent, and the conductivity (7^/3 can 
be calculated from the converged solutions for the self energy J]{uj,T). The dc conductivity 



is calculated as 12 



f —df 1 

fc,(T 

where Cxx = e^n/ {2hao),ao is the lattice spacing, / = [exp^u/T) + Va = 9e^/9fca, and 
A{k,u) = — ^ Im[uJ + fi — — S(c<j)]~^. The transverse conductivity is given by ||12| 

-df, 1 



axyiT) = BcxyJ du{^)^Y.^xy{k)A\k,u) (4) 



where c^y = 2|e|'^7r^ao/(3/i^), Bxiduxyik) = v1{dvy/dky)—VxVy{dvx/dky). The Hall coefficient 
Rh is given by axy/Bal^. 

We first discuss the "pinning" of the vHs in the total DOS with respect to variations 
in doping. One can parametrise the spectral function for — > 0,T ^ as A{k,uj) ~ 
r/((m*w + Ji- e^)2 + TT^r^) where T = -(l/7r)/mS(cj = 0), m* = 1 - {dJ:R/dtu)^=o 
and ft = fJ'{S,T) — S/j(5, T, u; = 0). The Luttinger sum rule ensures that fi{6,T = 0) = 
fi{S,T = 0,U = 0) = /Uo(^)) the non-interacting band chemical potential. At sufficiently low 
temperatures, T/D « 1, so for calculating the small u DOS A{k, u) can be approximated as 
S{m*u!+fi—ei^). Since the DOS in the interacting system, p(co'), is given by —{l/n)ImGii{Lj~^) 
eqn [2] reveals that for u ^ 0, p{uj) ^ po{m*uj + Ji). For T 0, ft po, which directly 
shows that if the non- interacting system had the vHs at the Fermi level {u = 0), then it 
survives at arbitrary U as long as m* is finite. The form of p{uj) also reveals that close to 
optimal doping the shift of the vHs from the Fermi level with doping is reduced by a factor 
of m* with respect to the non-interacting problem (c/. Fig. 1). As we will see later, it makes 
a significant difference to the transport that the vHs continues to be close to Ep for hole 
doping close to optimal doping, while for the electron doped [6 = —0.2) case the vHs in the 
DOS is ~ OMD below Ep. 
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The vHs in the cuprates is inferred from the saddle point observed in the ARPES " disper- 
sion curves" . The saddle point and the accompanying "flat bands" are ubiquitious across the 
cuprates 0, having been observed in the Bi 2212 [0 and 2201 fg, YBCO and NCCO 
|T^. However, comparison of LDA calculation and ARPES data |]I6|, along the symmetry 
lines r (0,0) ^ X , ^ d (^,0) ^ P (0,0), reveals that the measured "bandwidth", 
particularly along the X — Gi (zone corner - edge center) direction is reduced by a factor of 
2 with respect to the LDA result. Fig. 2 shows our results for the "dispersion curves", as 
inferred from the peaks in our A{k, to), for 5 = 0.2 (hole doped case) and 5 = —0.2 (electron 
doped case). For both cases there is a substantial reduction of the bandwidth, by factors 
of ~ 3 and ~ 4 respectively. In the optimally doped case, the "flat band" is at Ep and 
disperses less than ~ 10 meV along Gi — P. In the electron doped case the "fiat band" is 
~ 60 meV below Ep and disperses much more (along Gi — P). Also worth noting is the 
asymmetry in the narrowing of the dispersion relations for electron and hole excitations. 
For the hole doped case, the hole dispersion is narrowed more than the electron dispersion 
as the former correspond to the lower Hubbard band which contains the Fermi level. In the 
electron doped case, the structure is reversed, as the Fermi level lies in the upper Hubbard 
band 

The singular structure in the DOS close to Ep has non-trivial consequences for the T 
dependence of the damping rate P and the resistivity p^c- Although T,{lj, T) is "Fermi-liquid 
like" for T, — 0, with P and pdc ~ T"^-, beyond a low T scale that seems to be set by the 
proximity of the vHs to Ep the T dependence becomes "anomalous" as shown in Figs. 3 
and 4. In Fig. 3 we contrast the P calculated with the vHs to that for a symmetric, regular 
DOS, Po(e) = (2/7rZ})y^l — {e/Dy. For both the vHs and the regular DOS P starts out as 
T^, but in the vH case, and for hole doping it rapidly changes over to a quasi-linear form. 
In fact in the regime T ~ 0.01-D — 0.08D, at optimal doping P ~ T. The same effect is also 
seen in the T dependence of pdc, as shown in Fig. 4. We have not yet carried out enough 
systematic studies to pick out the dependence of the crossover scale on the parameters of 
the model. 

The weak coupling vHs scenario for anomalous transport has been criticised by Hlubina 
and Rice [|1^ on the basis that in this scenario 1/r is strongly k dependent. It is anomalous 
(1/rg ~ T,ej^ etc) only for k regions close to the vH saddle points (±^, 0), (0, ±^), and 
hence their contribution to p^c will be "shorted out" by the k regions away from these 
points, which have a normal, smaller 1/r (1/t^ ~ T^,e|). This criticism clearly does not 
apply to the above strong coupling results, where P is uniformly anomalous over the entire 
FS. While the T dependence of P and pdc are rather different from those for the regular 
DOS, the results for Rh are not dramatically different |TT|, except that now the (T = 0) 
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noninteracting Rh is also positive for the hole doped cuprates and a change in sign in Rh 
(at low doping, low T) occurs only for the electron doped cuprates. The local approximation 
does not reproduce the T dependence of Rh observed in the cuprates. 

In conclusion, we have shown that even in the context of a simple local approximation, 
the proximity of a vHs to the Fermi level leads to interesting strong correlation effects in 
ARPES and pdc that bear many similarities to the data in electron doped and in optimally 
or over (hole) doped cuprates. Recent ARPES data on under(hole) doped cuprates have 
been interpreted as showing evidence of a non-Luttinger FS, with regions of the FS near the 
vH points becoming "gapped" |T9|-^. There is no way one can access such features within 
a "local" approximation, with a k independent S. One can hope that putting back the k 
dependence of E, and of the vertex functions as "perturbative" , "non local" corrections to 
the above theory may cure the failings pointed above. 
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FIG. 1 The interacting DOS at low frequency. For the "optimally doped" {S = 0.2) case 
the vHs remains at Ep. For 80% doping variation the peaks are within a width of ~ 0.22. 
The corresponding width in the non interacting system is ~ 0.42. 

FIG. 2 "Dispersion" for complementary p and n doping; U = 4:,T = 0.01. The "flat 
band" is right at Ep for 6 = 0.2. The non interacting dispersion goes from —1 to 1 for the 
same scan. 

FIG. 3 Single particle damping rate for U = 4. Inset; P for the symmetric, regular DOS. 
FIG. 4 Dc resistivity, in units of e^7r/(2/iao), for representative doping. 
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